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Analytic Approach to the StabiUty of Standing Accretion Shocks: 
AppUcation to Core-Collapse Supernovae 
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^ ■ ABSTRACT 

We explore an analytic model of the accretion shock in the post bounce phase 
On , of a core-collapse supernova explosion. We find growing oscillations of the shock 

in the / = 1 and / = 2 modes, in agreement with a variety of existing numerical 
^ . simulations. For modest values of the ratio of the outer accretion shock to that of 

\0 . the inner boundary to the shocked flow, the instability appears to derive from the 

growth of trapped sound waves, whereas at higher values, postshock advection 



(N 



O , clearly plays a role. Thus the model described here may relate to the different 

Q . mechanisms of instability recently advocated by Blondin & Mezzacappa, and by 



9"' 
6 



Foglizzo and collaborators. 

Subject headings: supernovae: general — hydrodynamics — instabilities — shock waves 



• ■ 1. Introduction 

rS 

H , 

5r , There is a growing consensus that a large fraction of core-collapse supernovae, pos- 

sibly including the whole subset of Type Ib/c events, undergo significantly asymmetrical 
explosions. While clearly rotation and magnetic fields may play some role in this, instability 
associated with the stalled shock that forms at a radius of 100-300 km and lasts for a few hun- 
dred microseconds post bounce offers the intriguing possibility of generating asymmetrical 
explosions from a symmetrical progenitor, since such instability is usually dominated by low 
mode / = 1, 2 oscillations of the shock front. The original speculation for the origin of such 



^E. O Hulburt Center for Space Research, Naval Research Laboratory, Code 7674L, Washington DC 
20375 

lamingOnrl . navy . mil 



-2- 



osc illations (Herant et aljl992 : Herant 1995 ) in terms of co nvection behind the shock, follow- 



Chandrasekharl ( 1961 ) has been revised in recent years. Foghzzo. Scheck &: Janka ( 2006 ) 



mg 

demonstrate that post shock advection can stabihze convective instabihty (j Chandr asekhar 
( ll96ll ) only deals with a static case), and that convection alone cannot produce the dom- 



inant I = 1,2 modes. I n adiabatic simulations of gas with polytropic index 4/3, which as 
shown by iJankal (120011) is a good description of sh ocked inner regions of a core-collapse 



supernova, iBlondin. Mezzacappa fc DeMarind (12003 ) were the first to associate the insta- 
bility with the shock itself. iBlondin fc Mezzacappal ( 20061 ) attributed this to the growth 



of trapped sound waves, travelling essentially laterally around the central accretor. This 



must be a purely shock driven instable 
different to th e origin al sp eculations ( 



ity, since coiiyection is absent, an d hence qualitatively 



..Herant et al.lll992l : lHerantlll995l ). On the other hand, 

Foghzzo et al.l (120071 ) and iGalletti fc Foglizzo J2005h ~ attri bute this shock-driven instability 



to an "advective-acoustic" cycle, originally introduced by iFoglizzd (120021 ) in the context of 
black hole accretion flows. Vorticity perturbations generated at the distorted shock front 
advect inwards, and couple with outgoing sound waves at the inner boundary of the post- 
shock accretion flow. Upon reaching the outer shock, these sound waves further distort the 
shock front , leading to a p ositive feedback. While the mechanisms of IBlondin fc Mezzacappa 
J2006h and iFoglizzd (120021) appear to produce similar outcomes in non-rotating cases, indeed 



Foglizzo et a 



(120071) actually confirm the numerical values of frequencies and growth rates 



measured by IBlondin fc Mezzacappal (120061 ) . it is possible the presence of rotation might lead 
to substantial differences, which we discuss further below. 

Instabilities of the accretion shock also are observed in simulations with a more ac- 
curate treatment of supernova microphysics, with a realistic equation of state, neutrino 
transp ort and attendant n uclea r reactions, and mo re physically motivated boundary condi- 
tions. IScheck et al.l (120041 ) and IScheck et al.l (120061) investigate the role of such instabilities 



in generating pulsar natal kicks, while iKifonidis et al.l (l2005r) study t he ejection of metal 
clumps to the outer regions of supernova ejecta. iBurrows et al.l (120061 ) also see the stand- 
ing accretion shock instability 100 - 300 ms postbounce. The important feature of their 
simulation, though, is a core gravity wave at late times that transfers acoustic power to 
the shock to continue powering the explosion, the initial perturbati on for which is possi- 



bly seeded b y the accretio n shoc k instability. IBurrows et al.l (l2006l ) also favor the inter- 
pretation of iFoglizzo et al.l (120071 ) of an advective-acoustic cycle, though in this and the 
other more realistic simulations cite d above, the exact mechanis r n of i nstability is difficult 
to identify. To avoid this problem, lOhnishi. Kotake fc Yamadal (120061 ) performed numeri- 
cal simulations based on unper turbed spherically symmetrical shock accretion flows given 



bv lYamasaki fc Yamada J2005h. a s ort o f "mi ddle ground' between th e adia batic flows of 
Blondin. Mezzacappa &: DeMarind (120031 ) and IBlondin &: Mezzacappal (120061 ) and the full 
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supernova simulations of other workers. lOhnishi. Kotake fc Yamadal (120061 ) clearly identify 
and measure growth rates for / = 1 and 2 mode instabilities. The real parts of the frequencies 
scale with JJ^" (l/cg + 1/fr) dr, the integral between the inner radius of the accretion flow 
and the shock of the sum of the inverses of the sound sp eed and advectio n sp eed, suggesting 
that a n advective-acoustic mechanism, as advanced by Foglizzd (120021 ) and iFoglizzo et al. 



(120071 ) , is at work. However lOhnishi. Kotake fc Yamadal (120061 ) only perturb the radial ve 



locity component of the initially spherical symmetrical unperturbed flow, which may unduly 
influenc e their results, compared to a llowing perturbations in polar and azimuthal directions 
as well. lYamasaki &: Yamadal (|2006l ) include rather more physics, and again based on the 
frequencies at which modes grow, conclude that an advective-acoustic cycle is at work. 

In this wor k we p resent another approach to the problem. Following the work of 
Vishniac &: Ryul (119891 ) on planar shocks in the interstellar medium, we derive an approx- 
imate dispersion relation for oscillations of the accretion shock by adding the effect of a 
gravitational field and spherical geometry. In the case that the postshock advection is ne- 
glected, a quartic equation (actually a quadratic in uj"^ where u is the frequency of oscillation) 
results, becoming a quintic equation in uj with the inclusion of advection. We find that in 
both cases for 7 ~ 4/3, for all except the highest values of the ratio of the shock radius 
to the radius of the inner boundary, modes with / = 1 grow the fastest (in fact for cases 
where the postshock radial velocity, Ur = 0, no other nonradial modes grow). This sug- 
gests that for / = 1 in 7 = 4/3 gas in this regime, postshock advection is not crucial to 
th e operation of the instab i lity. Thus in general these results will support the conclusions 
of iBlondin fc Mezzacappa J2006h that the instability proceeds by the growth of trapped 
waves in this regime. Nevertheless, regions of / and 7 parameter space exist where we find 
modes that are only unstable in the presence of postshock advect i on, w hich we interpret as 



an ad vective-acoustic cycle, as suggested by iGalletti fc Foglizzd (120051 ) and iFoglizzo et al. 



(120071 ). In both cases, the frequencies of the growing modes are similar, suggesting that 
frequency alone is not a good discriminator of the mechanism of instability. As will be seen 
below, this might be expected in a comparison of a perturbation advected radially between 
the shock and inner boundary, and a sound wave traveling essentially laterally around the 
shock. 



2. Analytic Theory 

We consider an unperturbed model in which spherically accreting plasma is decelerated 
at a spherically symmetrical shock, before accreting onto a protoneutron star. The accre- 
tion shock is at radius r^, and we take an inner boundary at r^, maintained at constant 
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pressure, where plasma cools and decouples from the postshock flow. The postshock flow 
is modeled by the Bernoulli equation, as given in Appendix A. Gas with polytropic index 
7 < 1.5 decelerates, and gas with 7 > 1.5 accelerates away from the shock towards the inner 
boundary. We shall be solely concerned with 7 < 1.5. The dominant variations with radius 
are in pressure and density, as shown in Figure 1. 



To treat the perturbation, we follow in large part the methods and notation of IVishniac fc Ryu 



(1l989l ) who derive an approximate dispersion relation for application to shocks with arbi- 
trary postshock structure. In spherical coordinates, we write the hydrodynamic equations; 
an equation of continuity and momentum equations in the r, 9 and directions as 



Continuity: 

d{p + Ap) 
dt 

+ (P + Ap) 



,^ .5^P + Ap) t^9^p + Ap) (u^ + t;^)9(p + Ap) 



1 d sin 9v0 1 



dr r 89 r sin 9 d 

d{u^ + v^) 1 d {r"^ {Ur + Vr)) 



r sin 9 89 r sin 9 



dr 



Radial Force: 

d {Ur + Vr 



(P + Ap) 



dt 



d{Ur + Vr) f U4, + V4,\ d {Ur + Vr) {u^ + V^)" 
\ yUr ' Vr J r^ ~r 



dr 



r sin 9 



d{P + AP) d 
dr dr 



GM {p + Ap) 



(2) 



Poloidal Force: 

dv. 



(p + Ap) 



dVe VgUr „M<^ 

H cote*— 

dt dr r r 



-2cot^^ + 



dvg 



r r sin 9 dq 

ld{P + AP) 

r d9 



(3) 



Azimuthal Force: 



(P + Ap) 



d {u, + v,)_ ^ ^ d{u, + v,) ^v_^du^^ ^u,ve 



+ (P + Ap) 



dt dr 

+ V^ \ d (m0 + V^) {U^ + V<p) {Ur + Vr) 

r sin 9 



r d9 r 

1 d{P + AP) 
r sin 9 dd 



(4) 



where p and Ap are the initial and perturbed density respectively, Ur and Vr are the initial and 
perturbed radial velocities and and f ^ and the initial and perturbed azimuthal velocities. 
The initial poloidal velocity is zero, and its perturbed value is fg. P and AP are initial 
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and perturbed pressures, and GM is the gravitational constant times the enclosed mass 
which is dominated by that of the protoneutron star. Writing Ap/p = 5, and assuming 
5 oc exp {iojt + J \dr) Yi^ {9, (p) where A = A (r) we derive a linearized continuity equation, 

1 d sin 9v0 



(12 

lUJ + Vr \ — H 

L r 



imVfp dvr 
r sin 6 dr 



r sin 6 89 



0, 



(5) 



where uj' = uj + mu^fi/r sin 9 = uj + mQ where Q = u^/r sin 9 and L = (1 — 7) r is the density 
scale length, defined by din p/dr = 1/L. The derivatives of terms in and have been 
taken to yield the m dependence, but the 9 derivative of the term in vg remains explicit. 
Below (in section 3.2) we will assume oc rsin^ rendering Q constant. 



The radial force equation is linearized as follows. The zero order equation gives 

dUr dUr UA) dUr 

+ u 



p- 



dt 



+ 



dr r sin 9 



_dP_ d_ f GMp \ 
dr dr \ r J 



(6) 



If dur/d(j) = and dur/dt = 0, and 7 = 1 (ur = 0) or 7 = 1.5 {dur/dr = 0), then 
equation 6 evaluates to zero. We will assume that this is approximately true for all 1 < 
7 < 1.5. (Certainly Urdur/dr < throughout this range, whereas for 7 > 1.5, Urdur/dr > 
0). This helps ensure that for all cases where the postshock gas decelerate s, the radial 
/ = mode will be sta ble, in agreement with the conc lusions of iNakayamal (Il992l ). and 
Burrows fc Goshyl (119931 ) and lYamasaki fc Yamadal (120051 ). where we are assuming a neutrino 
luminosity below the critical value in these last two references. Taking AP = P6 (i.e. an 
approximately isothermal perturbation; this gives a considerable simplification, see Appendix 
B) the equation to first order in small quantities becomes 



d_ 

dr 



'Xlirri V f 



VQ dUr 

r de 



+ 6 



exp J i^dr 



2ui 



+ 



exp J i^dr. 



(7) 



We note that 



d 



1 



dr \_iuj'/ur + A 



exp / A + iu' / Urdr 



—dX/dr + i (uj/u^.) dUr/dr 

{iuj'/Ur + A)^ 



+ 1 



exp / X + iuj' /urdr. 



For a global mode we expect A to be of order 1/r, (in fact for 7 = 4/3, A ~ 2/r; see equations 
16 below). If UJ represents laterally propagating sound waves, u is in the range Cs/r — Cg/ \L\, 
whereas u 2tt {r/ur + r /cg) for radially propagating vorticity perturbations. In either case, 
with dur/dr = (ur/r) (3 — 27) / (7 — 1) from Appendix A 
and we may put 



dA/dr+^(a;/^^r)(3-27)/(7-l) . 

{iuj'/Ur + Xf 



d 

dr 



1 

iu'/Ur + X 



exp / A + iuj' /urdr 



exp / A + iuj' /urdr. 



(9) 
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In equation 7 we neglect terms in u^/cl, u^/cl (of order (7 — 1) /2^) and integrate using 
equation 9, further assuming \iu!'/ur + A| >> |l/r|, to get 

2u^v^ ve dur , S [ GM ^\ 

Vr = -T^T^ - —, 7^ + —, tA C. . (10) 

luj'r loj'r oO loo' + UrX \ r J 
With = and dur/d9 = 0, equation 7 gives 

—Ur = I — Cj A - lUJVr - -^Vr = I "^'^ ~ ) ^^^^ 

where we have substituted from equation 10 for Vr to simphfy the right hand side. Taking 
Ur oc r^^'^'^y^'y-'^) from Appendix A, we derive 

Expressions for vg and are easily derived in the limits — > or ^ 0. Staying with 
the nonrotating case for the time being, we give for u^j) — 0: 

" r do iui+Xur 

y ^_,JB£L (13) 

where the primes have been dropped from uj. 

Substituting our expression for v^, dvr/dr, vg and v^p into the linearized continuity 
equation (5) we find 

. , 6X [GM 2WI 47-51 A cl I {I + 1)6 ^ ,,,, 

luS + — T + — — + A + UrXS + ^- ^ — = 0, (14 

luj + AUr \ r J \^ 7 — Ir / lu -\- XUr 

where the m dependence disappears (as it should in spherical symmetry) using properties of 
the spherical harmonics. This can be solved for A to give 

~ \2L ~^ 7-1 2r J GM/r-cj 

I /j_ I 47-5 /-I I ^ iUrUJ I A Uj'^-l{l+l)cl/f 

'^\2L^ 7-1 Y ^ ^(GM/r-c2)(l/L+(47-5)/(7-l)/r-) ^ ^ (l/L+(47-5)/(7-l)/r-)"(GM/r-c2) " ^^^""^ 

We write (neglecting terms in u^/cl) 



GM/r-c 



.2 



K + X- =-(i + ?Efi)-ai^. (16) 
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As discussed by Vishniac fc Ryul ( 1989 ). a third solution exists, with 5 = everywhere. 
In that case Vr oc exp (— z J u/urdr) /ur, and vg and oc exp (— i J u/urdr) /urr. The 
hnearized continuity equation (5) then gives 



1 2 
L r 



or 



0. 



(17) 



Taking dvr/dr = —iVrUj/ur » Vr {^/L + 2/r) we find 

Vr = ■ V±. 



The derivation of the dispersion relation from the boundary conditions in terms of 
these three solutions is carried out in Appendix B, and given by equation B6. We use 
= A-t/ (1 + X±u/iuj), with A^ and Aj evaluated at the outer shock and inner boundary 
respectively in terms of Urs and Uri, the advection velocities at these locations (see Appendix 
B for more details). Then in equation B6 we approximate 



j_ _i_ 47-5 1 

2L 7-1 2r 



Qj 



2L 



+ 



47- 



' GM/r 
-5 1 



1 i+fsQ 



ui GM/r-cj 



-1 2r 



1 + Qf^) + 



1 



+ 



2uj y2L 

im)cVr^ 

GM/r-cl Q l-f3Q uj GM/r-cl 



47-5 1 
7-1 2r 



47-5 1 
7-1 2r 



(19) 



and 



l{l + l)cl/T^-u'^ 

GM/r-cl, 



1 I 47-5 1 

L 7-1 r ) 2(GM/r-c|) 



47-5 1 
7— 1 r 



2lu 



1 + Qti^ 



2uj 



Ur 



2 - Qizfcj + ^13^) (2 + Qjz^ - ^T^) (20) 



1 1-13Q 



where we have put 

Q' = Jl + 4i 



(GM/r-c2)(l/L+(47-5)/(7-l)/r) 



+ 4 



t^2-;(f+l)e2/r2 



(l/L+(47-5)/(7-l)/r)"(GM/r-ci) 



(GM/r-c2)/(l/L+(47-5)/(7-l)/r) 



and neglect terms of m^/c^ and higher. These then give 



(21) 



u;^ + u;^[^Urs{i + i)-^-Yr]VT^ 



1(1+1) 



GM/r-cl ( 1 
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i ( GM 
'2 \ r 



I) 



1 _|_ 47-5 1 
L 7-1 r 



{-rs {{ + ^) (1 - Q'^) + ^ (1 



+ : 



GM/r-cl 
aL 

GM/r-cl 
2al 



2aL 



Y3g^ -I- fU-ri -I- Q ) 



GM/r-cl 
2aL 



1 _|_ 47-5 1 
L 7-1 r 



i _|_ 47-5 1 

L ~'~ 7-1 r 



— 



?(< + l)-»sMrs 



L 7-1 r 



a + A.) 

KL ^ 2rJ 



lUr 



1 + Qt 



+ i 



aL 



(22) 



If Ur = 0, equation 22 is very similar to equation 13a of IVishniac fc Ryul (Il989l ). and is 
identical if we take r — 00, L —L (to agree with their sign convention), / (/ + 1) /r^ — A;^ 
and 7^1 (i.e. VgUr — > cj.). In their case the l imit Q ^ 1 giv es a dispersion relation similar 



to the exact solution in the thin shock limit (jVishniad 119831 ). and identical to this limit if 
terms describing the evolution of the shock are neglected from the exact solution. They then 
proceed to argue that taking Q 7^ 1 gives a suitable approximate dispersion relation away 
from the thin shock limit. Our case is pa rticularly suitable since t he standing accretion shock 
can be treated as an equilibrium state (iBurrows fc Goshylll993l ). i.e. it does not evolve in 
our approximation. 

Just t aking Ur = 0, equation 2 yields a quadratic equation in u;^ for which simple solu- 
tions exist. IVishniac fc Ryul (119891 ) find growing shock oscillations in the case of decelerating 



shocks, with a post shock rarefaction. In our case, the growing oscillations exist for density 
increasing with distance behind the shock. This difference arises because the square of the 



sound speed, cj. 



GM/r in our case, and GM/r > c^. When Ur = (except in terms 



l{l + l)VsUr 

2r2 



( GM _ 2\ (J_ I 5_\ 47-5 1 
\ r sj \2L^ Ar) 7-I r 



in VgUr) and taking Q \ 

ijj 

( GM _ 2\ ]_ (]_ I 5_\ lil+l)v8Ur 
\ r ^s) L \L~^ 2r) r'^ 



where 



W 



l{l + l)VaUr 



When 1 >> W, this gives a pair of sound waves for given / with 



/ (/ + 1) VsUr/r^ - {GM/r'^ - c^/r) (1/L + 5/2r) 



47-5 
7-1 



taking the positive sign, and the negative sign gives a pair of gravity waves with 

( GM _ 2\ ]_ ( ]_ I 5_\ l{l+l)vsUr 
\ r ^s) L\L^ 2r) r'^ 



i(/ + l)fsMr 



( GM _ 2\ (I I 5_\ 47-5 1 1 
V r ^sj \l 2rJ 7-I rj 



(23) 
(24) 

(25) 
(26) 
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where a is a constant coming into our definition of the interior boundary condition, as 
explained in Appendix B. For 1 < W, we have four roots of the form uj = ± {f + ig) and 
uj = ±{f — ig), and all roots have a combined character of acoustic and gravity waves. 
Reinstating the terms in odd powers of u with nonzero Ur, the roots separate into pairs 
of what may be termed acoustic and gravity waves, although the difference in frequency 
between them is generally rather small and both pairs are of presumably mixed character, 
and one imaginary (i.e. purely growing or damping) mode. Assuming a value of Q, we solve 
the quintic equation numerically. Results for each growing mode are then iterated until the 
value of Q produced by the frequency and wavenumber matches the input. Results of this 
procedure for growing modes are given in Table 1 for 7 = 4/3, a = 1.5, and in Table 2 for 
7 = 1.36, a = 2, the values of a being chosen to reproduce as far as possible the stability 
properties of the radial mode. We compare in each case the results of the full quintic equation 
with a simpler case taking Urs = Uri = 0. 

With 7 = 4/3, taking GM/rcl ~ 1/ (7 — 1) = 3 (from equation A2 neglecting u^), and 
without advection, we find growth for / = 1 for all values of Ts/vi < 0.1. Growth at Z = 1 
at r^/rj —i> cx) and at / = 2 for all r^/rj require non-zero advection terms in the dispersion 
relation, though the / = 2 mode is stable as r^/rj — >■ 00, even in the presence of advection. 
The / = mode is marginally stable in the absence of advection. With advection included, 
it is stable for the lowest r.^ /r,:, but otherwise unstable. T aken together, there is reasonable 
quahtative agreement with iBlondin fc Mezzacappa] (120061 ) . but quantitatively, the model is 
less stable than it should be. Growth rates are too high, and the / = mode becomes 
unstable too quickly as r^/rj increases. Better qualitative and quantitative agreement is 
found in the model with 7 = 1.36 and a = 2 given in Table 2., and also illustrated in Figure 
2. The stability of / = is improved, and advection becomes relatively more necessary to the 
instability of I = 1. These results have asimple i nterpretation in terrn s of the mechanisms 
put forward by lBlondin fc Mezzacappal (120061 ) and lFoglizzo et al.l (120071 ). Where advection is 
not necessary to the instability, i.e. for the / = 1 mode at low values of r.^/r^, the wave growth 
must be due to trapped sound waves similarly to the simulations of iBlondin fc Mezzacappa 
( 120061 ). The inclusion of advection in these cases generally increases the growth rate, either 
due to a separate contribution of an "advective-acoustic" nature or because of the effect of 
postshock advection on the propagation of sound waves. In cases where no instability occurs 
unless ad vection terms a r e pres ent, the growth must be due to an advective-acoustic cycle 
similar to iFoglizzo et al.l (120071 ). In our model, these cases include all the / = 2 modes, as 
well as the / = 1 modes for large values of r^/rj. 



These results and interpretation add a new twist to the conclusions of iFoglizzo et al. 



( I2OO7I ). These authors demonstrate that for r^/rj > 10 the advective-acoustic cycle is the 
cause of the / = 1 instability. Their analysis makes use of a WKB approximation that 
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is not valid for lower r<j/rj. However they proceed to argue that since growth rates and 
frequencies vary smoothly as one goes to lower values of r^/rj, the mechanism of instability 
should remain the same. However our results suggests that even though the eigenvalues vary 
smoothly, the mechanism of in stability changes to more closely resemble that suggested by 



Blondin fc Mezzacappal (120061 ) . 



3. Discussion 



3.1. Preamble 



A number of simplifying assumptions have been introduced into the analysis to make 
the algebra tractable. We have assumed an isothermal perturbation, used a greatly sim- 
plified interior boundary condition, and have dropped all terms in the advection velocity 
of order higher than Ur/cg. Ne vertheless our model ap p ears t o capture the basic features 
of the numerical sirn ulati ons of iBlondin fc Mezzacappal (120061 ). and the analytic theory of 
Foglizzo et al.l (120071 ) and lGalletti fc Foglizzd (120051 ). in the appropriate regions of parameter 
space. Looking at the values of A± given by equation 15, we find when tu^ ~ / (/ + 1) c^/r^, 
|A_| >> |A+|, A_ > and A+^O for Q^l- Consequently, A_ gives a density perturbation 
5 = 5p/ p that is greatest at the outer shock and considerably smaller at the inner boundary, 
while A+ gives a density perturbation that varies much less with radius, but can be largest 
at the outer shock or at the inner boundary. When / = 0, A+ < and |A_| ~ |A+|. In this 
case the radial mode can have a large maximum of 5 at the inner boundary. Given the ad 
hoc nature of our inner boundary condition (see Appendix B), we might expect less accurate 
results for the radial mode than for nonradial modes, which does indeed seem to be the case. 
The vorticity (calculated from equations 13) increases inwards approximately as The 
third solution with 5 = everywhere must represent a purely vortical perturbation. 



3.2. The Mechanism of Instability 

So far we have developed a dispersion relation from the equations of hydrodynamics, and 
found instabilities through the existence of complex roots. We have not made any comment 
on the precise mechanism of the instability, other than to comment that sound /gra vity 
waves appear to grow. The vortical-acoustic cycle has been explained by iFoglizzd (120021 ). A 
distortion of the accretion shock produces a vorticity perturbation that is advected inwards. 
Upon reaching the inner boundary, the vorticity perturbation produces an upward moving 
sound wave, which reinforces the shock distortion when it arrives there, producing a net 
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increase in the original perturbation. Over many cycles, exponential growth ensues. 



Sound waves produced at the outer shock, in a realistic model, will not propagate 
towards the center due to the ef fects of refraction. In s tead, they will propagate around 



the circumference of the shock. iBlondin fc Mezzacappal (120061 ) speculate that sound waves 
produced at one position on the shock surface by a density inhomogeneity will propagate 
around until they meet on the opposite side, where their excess pressure produces a shock 
distortion that sends another pair of sound waves back again. It is not immediately clear 
without detailed calculation that such a process should produce a growing (as opposed to 
damping) oscillation. We speculate that growth must also be aided by the flow of plasma 
in through the accretion shock. In Figure 3 we show fiowlines for shocked plasma in an 
accretion shock with a 10% / = 1 modulation of the shock radius. When the shock is shifted 
up by the instability, incoming plasma in the equatorial regions is diverted downwards by the 
now oblique shock. This non-radial flow of shocked plasma enhances the oscillation and leads 
to furt her wave growth. This aspect is simil ar to the th in shell overstabi l ity of decelerating 
shocks (IVishniadll983l : IVishniac fc Ryulll989l ). (although I Velikovich et al.l (120051 ) put forward 
a different view), though in these cases the pressure produced by gravitational confinement 
plays no role, and overstability exists for a postshock density gradient oppositely directed to 
our case. Thus both the advective-acoustic (or vortical-acoustic) instability and the growth 
due to trapped sound waves have their origins in perturbations of either vorticity or pressure 
in the post shock flow produced by the distorted shock front. 



3.3. The Effect of Rotation 



There is considerable interest in the effects of rotation on core-collapse explosions. In 
particular, in which direction with respect to the rotation axis should one expect the strongest 
growth? This would have obvious implications for the directions of pulsar kicks, assuming the 
such instabilities are the correct mechanism. In the case of an advective-acoustic instability, 
we should expect the decrease in the postshock advection in the plane of rotation to give 
stronger instability and hence strongest growth along the axis of rotation. 

In the case of instability driven by trapped sound waves we can sketch out the effect of 
including a nonzero u^p in the forgoing analysis. We derive the extra terms to be included 
in the dispersion relation in the hmits — > and u'^/c^ — > 0. We also neglect any distor- 



tion of the shock front by the rotation. Such considerat ions require more detai 



neutrino cooling near the inner boundary at a minimum (IBlondin &: Mezzacappa 



ideally also neutrino lumino sities and mass accretion rates (see e.g. iBurrows fc Goshylll993 



ed ph ysics, 
2006h . and 



Yamasaki fc Yamadal l2005l ). and are beyond the scope of this paper. This keeps uq = 0, 
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which should be true in any case if Ur = 0, but even in the absence of advection, a distor- 
tion of the shock front may change the boundary condition leading to equation B2. The 
perturbed velocities are 



5X_ (GM_ _ _ 2m.Q,cl8 



mc^s^ I 2n5X sine ( GAI _ 2^ _ 2nc^ cos d d& 



With 



dvr 6 {GM/r - cl) X' 

— \ — ; , y^o) 

or r luj' 



the linearized continuity equation becomes 

after dividing through by 5, with solutions for Aj- 



A+ + A- = -(i + i) + 



2Qti 



^+^- GM/r-c'i rLu' GM/r-cj ' y'^^) 



We again have Vr oc exp j —iuj'/urdr if 5 = 0, so equation (14) still holds. Then, with 
unchanged boundary conditions, the dispersion relation is the same as before (with Ur = 
and u ^ u' in equation B6), but with the extra terms (neglecting the change 1/L + 
(47 — 5)/(7 — l)/r^l/L + l/r, which in any case is trivial for 7 = 4/3) 



l{l+l)vsUr mn ( GM /n l+/3^ A ( GM cl\ ( \ , 5 2mQ.cl /01 N 

r2 uj'r \ rL ^ ) \ ) \ tL L ) 2r ) rLu>> ' 

The new quintic equation is again solved numerically, with results given in Table 3 for the 
m = 1 component of the / = 1 mode. We generally find higher growth rates increasing with 
rotation rate for the m = 1 (m = remains unchanged), which would indicate higher growth 
in the plane of rotation. This happens even in cases where in the absence of rotation, the 
advective- acoustic instability seems to be dominant. We speculate that rotation increases 
the growth rate of trapped sound waves for nonzero m, so much so that they can grow even 
in cases with zero growth in the absence of rotation. 

At higher rotation rates, the distortion of the shock front cannot be neglected. We expect 
the equilibrium s hock position to be ob l ate, i. e. bulging out at the equator as demonstrated 



quantitatively by lYamasaki fc Yamadal (120051 ). We emphasize again that our results are not 



applicable in this regime. 
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Conclusions 



Using heuristic arguments based on IVishniac &: Ryul (119891 ) . and a few other approxi- 
mations, we have derived an approximate dispersion relations for oscillations of a standing 
spherical accretion shock. The pursuit of an analytic treatment further through the problem 
than has previously been attempted arguably allows greater physical insight, albeit at the ex- 
pense of reduced quantitative accuracy. While our model might not be considered much more 
than a "toy model" , it does capture most of the features of previous works of the instability 
of such shocks. It offers a poten tial r esolution to the mild disagreement between the works 
of iBlondin fc Mezzacappal (120061 ) and lFoglizzo et al.l (120071 ). We find results consistent with 
both works, but in different regions of parameter space. For low values of r^/rj, the growth 
of trapped waves dominates in the / = 1, and appears to be the most important mechanism. 
The inclusion of post shock advection gives a small change, usually an increase, in the growth 
rate. At higher I, or Vg/ri, postshock advection appears to become crucial to the instability, 
which leads us to interpre t it as an advective-acousti c cycle. In a treatment of 7 = 4/3 gas in 
core collapse supernovae, iGalletti fc Foglizzd (120051 ) only find adv ective-acoustic instab ility 
for Ts/vi > 3.5, with the growth rate increasing out to r^/rj ~ 10. iFoglizzo et al.l (120071 ) are 
only able to demonstrate the operation of the advective-acoustic cycle for r^/rj > 10, due 
to the nature of the approximations involved, but conjecture that it should also operate at 
lower r^/rj. Our models agree with their demonstrations at large shock radii, but suggest 
that at smaller shock radius, the instability should change to resemble that observed by 
Blondin fc Mezzacappal (120061 ). 



The precise mechanism of instability probably becomes more significant in the rotating 
case. The advective-acoustic instability must grow preferentially along the rotation axis, 
where the postshock advection is not affected by the rotation. By contrast, we conjecture 
in section 3.3 that rotation should enhance the growth of trapped sound/gravity waves 
principally in the plane of rotation, leading to stronger growth in directions perpendicular 
to the rotation axis. Such effects may have important consequences for pulsar natal kicks, if 
these derive from hydrodynamic kick mechanisms. 



This work has been supported by the Chandra GI Program, the NASA LTSA and 
APRA programs, and by basic research funds of the Office of Naval Research. I acknowledge 
a scientific discussion with Ethan Vishniac, and an extremely careful and constructive reading 
of the paper by an anonymous referee. 
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A. Bernoulli Equation for the Postshock Flow 

Here we summarize the mode l for t he postshock flow in spherical accretion used by 



Blondin. Mezzacappa &: DeMarind (120031 ) . The flow is given by the Bernoulli equation 



2 27 P 2GM , 

ul + '— = Al 

7 — 1 p r 

and just behind the shock itself the jump conditions give the velocity, density and pressure 
as 



p _ / 2GM 1 /A9^ 



«^ = _7_1 /2GM 

^ 7±1 M r^rT- 1 
y 7-1 47r V 2GA/ r| 

2 M 

7+1 47r ' 

Evaluating P/ = constant from the jump conditions and substituting gives 

^ 47 / 7-1 \ " ^,3-27^/1-7 _ 1 = (A3) 



(7 + 1) (7-1) V7 + I 

where r' = r/rg and u' = Ur\/rs/2GM . Neglecting the flrst term (i.e the kinetic energy), 
Ur oc r^^~'^"'^^^'^^^^ and hence p oc r~^^^'^^^\ These results are exact for 7 = 5/3. In this 
approximation, the density decreases monatonically with increasing radial distance, and L 
deflned above is L = r (1 — 7). We note that for 7^1, close to the shock, the density may 
increase with increasing radial distance, i.e. a postshock rarefaction exists. For 7 = 1, this 
region extends between O.TSr^ < r < r^, and disappears for 7 > 1.1. 



B. Boundary Conditions and Dispersion Relation 



We recapitulate and modify very slightly the approach of IVishniac fc Ryul (Il989l ). Let 
B± be the amplitudes of the X± oscillations of the accretion shock. Then at the accretion 
shock, 6 (0) = —5rs/L — 56rs/2rs — 25vs/vs ^ —5rs/L — 55rs/2rs = — (1/ L + 5/2r) Vr/iuJ, 

which comes from the radial variation of the shock ram pressure. We follow IVishniac fc Rvul 

(11989 ) and neglect the term in Svs, which renders the perturbation isothermal. iBlondin fc Mezzacappa 



( I2OO6I ) also comment that the term in Svg is smaller by a factor of a few than the term in 



Srs in determining the shock ram pressure. Hence 



f{GM/r 



{B+X'^, + B_\'_,) + 



5 

2r 



(Bl) 
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where A'^ = \±/ (1 + \±u/iuj), and the subscript s indicates that this is evaluated at the 
outer shock. 

Also at the accretion shock, (r^) = VsV±_5rs = {vs/iuj)'V ±Vr {rg), which after some 
rearrangement gives 



Vs {GM/r - cfj , 



. = - ---- "-V i + B_,\'_) - ^Vi ( ■ t;! ) . 



:b21 



Eliminating V_l ■ v±_ between equations Bl and B2 gives 

(u'L' ( 1 - ^^^^) - - c^) AV,) + 



B_ uo'L' 1 







(B3) 



where V^-B± 
l/L + 5/2r. 



-/ (/ + 1) B^/r\ Vi V± ■ = -/ (/ + 1) • lil /r2 and l/L' 



The interior boundary, and hence the boundary conditions here, are less well defined 
in our model. We are assuming that the region where strong neutrino cooling sets in and 
causes the accreting mat ter to decouple from the hot postshock fiow will give a boundary from 
which waves can refiect. IVishniac fc Ryul (119891 ). considering a plane parallel shock, assume 
a constant pressure boundary condition, 5 = —5ri/L = —Vr/iujL. We assume something 
similar; 



5 (rj) = 5+ exp / A+dr + B exp / X^dr 



(r,) 6 in) {GM/r -cl) 



itoaL 



(B4) 



We do not take the local value of L in this boundary condition, but take the same value as 
assumed above at the forward shock, modified by a constant a which we vary to find the best 
match to the / = stability. In any case, 5 {ji) = — f,. (rj) /iuaL (r^) << —Vr (r^) /iuL (rj) 
as it must be, because the background density increases significantly over that predicted by 
the Bernoulli model due to radiative cooling by neutrino emission. We prefer to cast the 
interior boundary condition in terms of L and r evaluated at the outer shock, because of 
the algebraic simphfication it produces, and the somewhat easier con ditions requ ired t o kee p 



the / = mode st able. It is clear from the work o f var ious a uthors iNakayamal (e.g. 119921 ): 
Burrows fc Goshyl (e.g. Il993l ) ; lYamasaki &: Yamadal (e.g. l2005l ) that this must be true. These 
last two references, in particular, include considerably more physics, and are able to predict 
the shock radius, rather than just specifying it as we do, but they only treat radial stability. 
Consequently we have 



5+ exp / (A+-A_)c;r<^ 1 



GM 



+i 



-B^ { 1 



GM 



B5) 
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where the X±i are evaluated with u at the inner boundary, but otherwise with L and r at 
the outer shock. 

Now ehminating / B_ between equations B3 and B5 gives the dispersion relation: 



+ - c?)' (i + [K^'-s - AV.A'_,/5«] = (B6) 

where 

f3^ = exp r - (A+ - A_) dr = exp H + ^^L^-] Qdr ^ ( -] . (B7) 



L •y — 1 r J \r 
For 1 < 7 < 1.5, /3 — >• as Vi/rg — >■ 0, and as 7 — > 1 for fixed r^/r 



REFERENCES 

Blondin, J. M., Mezzacappa, A., & DeMarino, C. 2003, ApJ, 584, 971 
Blondin, J. M., & Mezzacappa, A. 2006, ApJ, 642, 401 

Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, ApJ, 640, 878 
Burrows, A., & Goshy, J. 1993, ApJ, 416, L75 

Chandrasekhar, S. 1961, Hydrodynamic and Hydromagnetic Stability (New York: Dover) 
Foglizzo, T. 2002, A&A, 392, 353 

Foglizzo, T., Scheck, L., & Janka, H.-Th. 2006, ApJ, 652, 1436 

Foglizzo, T., Galletti, P., Scheck, L., & Janka, H.-Th. 2007, ApJ, 654, 1006 

Galletti, P., & Foglizzo, T. 2005, in proceedings of the SF2A meeting, 17-30 June 2005, 
Strasbourg, EDPS Conference Series in Astronomy and Astrophysics, ed. F. Casoli, 



T. Contini, H. M. Hameury & L. Pagani, astro-ph/0509635| 
Herant, M., Benz, W., & Colgate, S. 1992, ApJ, 395, 642 
Herant, M. 1995, Physics Reports, 256, 117 



Janka, H.-Th. 2001, A&A, 368, 527 

Kifonidis, K., Plewa, T., Scheck, L., Janka, H.-Th., & Miiller, E. 2005, A&A, 453, 661 

Nakayama, K. 1992, MNRAS, 259, 259 

Ohnishi, N., Kotake, K., & Yamada S. 2006, ApJ, 641, 1018 

Scheck, L., Plewa, T., Janka, H.-Th., Kifonidis, K., & Miiller, E. 2004, Phys. Rev. Lett., 92, 
011103-1 

Scheck, L., Kifonidis, K., Janka, H.-Th., & Miiller, E. 2006, A&A, 457, 963 

Velikovich, A. L., Zalesak, S. T., Metzler, N. & Wouchuk, J. G. 2005, Phys. Rev. E, 72, 
046306 

Vishniac, E. T., & Ryu, D. 1989, ApJ, 339, 917 

Vishniac, E. T. 1983, ApJ, 274, 152 

Yamasaki, T., & Yamada, S. 2005, ApJ, 623, 1000 

Yamasaki, T., & Yamada, S. 2006, ApJ, submitted, astro- ph/060658l| 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-18- 



Table 1: Eigenvalues of growing modes < / < 2: 7 = 4/3, a = 1.5 





no 


advection 




advection 












I = 1 


Z = 2 


2 


±0.569 


±0.525 


±0.569 


±0.546 


±0.771 






-0.242i 


+0.034i 


-0.284i 


-0.326i 


3 


±0.391 


±0.443 


±0.447 


±0.520 


±0.730 






-0.253i 


-0.090i 


-0.305i 


-0.312i 


5 


±0.257 


±0.387 


±0.453 


±0.511 


±0.707 






-0.248i 


-0.155i 


-OMOi 


-0.2m 


10 


±0.145 


±0.342 


±0.479 


±0.516 


±0.693 






-0.238i 


-0.178i 


-0.311i 


-0.287i 


00 






±0.511 


±0.532 










-0.199i 


-0.311i 





Note. — Real and imaginary frequencies are in units o{\cs/L\. Growing modes have a negative imaginary 
part. 
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Table 2: Eige nvalues of growing modes 0</<2: 7 = 1.36, a = 2 





no 
Z = 


advection 


1^0 


advection 
I = 1 


Z = 2 


2 


±0.506 


±0.495 
-0.193i 


±0.523 
+0.010i 


±0.490 
-0.222i 


±0.829 
-0.173i 


3 


±0.347 


±0.404 
-0.175i 


±0.353 
±0.500i 


±0.454 
-0.228i 


±0.818 
-0.171i 


5 


±0.219 




±0.264 
-0.058i 


±0.425 
-0.220i 


±0.811 
-0.176i 


10 


±0.044 




±0.287 

-o.m 


±0.401 
-0.201i 


±0.806 
-0.194i 


oo 






±0.343 
-0.147i 


±0.383 
-0.119i 





Note. — Real and imaginary frequencies are in units of \cs/L\. Growing modes have a negative imaginary 
part. 
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Table 3: Eigenvalues oil — 1, m — 1^ — 1.36, a — 2 



rjr, 


0.0 


0.05 


U4, 

0.1 


0.2 


0.3 


0.1 


2 


-0.490 
-0.222i 


-0.498 
-0.231i 


-0.502 
-0.242i 


-0.511 
-0.264i 


-0.520 
-0.282i 


-0.527 
-0.298i 




0.490 
-0.222i 


0.489 
-0.205i 


0.484 
-0.190i 


0.473 
-0.154^ 


0.461 
-0.091i 




10 


0.401 
-0.201i 


0.430 
-0.224i 


0.433 
-0.241i 


0.441 
-0.264i 


0.449 
-0.281i 


0.456 
-0.294i 




-0.401 
-0.201i 


-0.425 
-0.171i 


-0.430 
-O.llSi 


-0.464 
-0.045i 


-0.483 
-0.032i 


-0.495 
-0.027i 



Note. — Real {lj' = u + mO) and imaginary frequencies are in units of \cs/L\. Growing modes have a 
negative imaginary part. For m = — 1, real parts take opposite signs. 
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Fig. 1. — Variation behind the accretion shock of the pressure, density, inward flow velocity, 
sound speed and gravitational potential for a 7 = 4/3 model. All quantities are normalized 
to unity at the accretion shock and radius unity, except the inward flow velocity, which is 1/7, 
the shock velocity being unity at radius —1. The pressure and density exhibit much bigger 
variations with distance than any other quantity, motivating us to simplify the problem 
by taking the sound speed, inward flow velocity and gravitational potential to be constant 
throughout the shocked plasma. 
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Fig. 2. — Angular frequencies and growth rates of unstable modes plotted against r^/rj. Solid 
lines give present re sults, (black without ad v ection ; grey with advection), and the dashed 
lines are those from iBlondin fc Mezzacappal (l2006l ). all in units of Cs/L evaluated close to 
the accretion shock. The growth rate of the / = mode without advection is 0; it is only 
marginally stable, and is not plotted. 
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Fig. 3. — Post shock advection velocity flow lines for an accretion shock with 10% / = 1 
modulation of radius for 7 = 4/3 (upper panel). The shock perturbed upwards (solid line) 
from its equilibrium position (dashed line) has its surface distorted such that incoming 
plasma flows preferentially towards the lower density and pressure region at the bottom of 
the figure. This helps induce the growing oscillation. 



